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ABSTRACT 



We have recently identified a new radial migration mechanism resulting from the overlap of spiral and bar resonances in galactic 
disks. Here we confirm the efficiency of this mechanism in fully self-consistent, Tree-SPH simulations, as well as high-resolution 
pure N-body simulations. In all barred cases we clearly identify the effect of spiral-bar resonance overlap by measuring a bimodality 
in the changes of angular momentum in the disk, AL, whose maxima are near the bar's corotation and outer Lindblad resonance. This 
contrasts with the smooth distribution of AL for a simulation with no stable bar present, where strong radial migration is induced by 
multiple spirals. The presence of a disk gaseous component appears to increase the rate of angular momentum exchange by about 
20%. The efficiency of this mechanism is such that galactic stellar disks can extend to over 10 scale-lengths within 1-3 Gyr in both 
Milky Way size and low-mass galaxies (circular velocity ~ 100 km/s). We also show that metallicity gradients can flatten in less than 
1 Gyr rendering mixing in barred galaxies an order of magnitude more efficient than previously thought. 
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1. Introduction 

In the past few decades, discrepancies in the solar neighborhood 
age-metallicity relation have implied that effective radial migra- 
tion (i.e., redistribution of an gular momentum^ must be taking 
place i n the MiUcy Way disk (Edvardsson e t al. 1993; Haywood' 
2008t ISchonrichandBinnevI 120091: see Minchev & Famaevi 



20101 for a comprehensive discussion). In parallel, there is now 



considerable observational evidence that the non-axisymmetry 
of the Galactic potential can cause significant perturbations in 
the motion of stars of all ages. Evidence for this c omes from, e.^ 



the moving groups in the solar neighb ourhood ([Dehne 
ILJ ' containing stars of very different ages dFamaev et al.L l2005h or 
. !^ ] the non-zero value of the C and K Oort constants for red gi- 
• ant stars in t he ex tended local disk (Oiling & Dehnen, 2003i 
$-H ' ISiebert et al.L l2010h . All of these can be expla ined by the ef- 
, feet o f reson ances associated with a centr al bar d Minch ev et aL 
2010[ l2007h or spiral s tructure (SS) (lOuille n and Minchev 
2005trOuillen et al.[|20Tol) . 



Until recently it was accepted that efficient radial mixing 
of st ars in galactic disks was cau sed solely by transient spi- 
rals (ISel lwood and Bi nnevl l2002l hereafter SB02). However, 
lOuillen e t al. (2009) showed that small satellites on radial, in- 
plane orbits can cause mixing in the outer disk and thus ac- 
count for the frac tion of low-meta lhcity stars present in the so- 
lar neighborhood (lHavwoodll2008h More over, we have recently 
demonstrated (iMinchev & Famaevi l2010i hereafter MFIO) that 
a strong exchange of angular momentum occurs when a stel- 
lar disk is perturbed by a central bar and SS simultaneously: 
our test-particle simulations allowed us to attribute this effect 



to the overlap or first and second order resonances of each per- 
turber. The presence of multiple patterns in g alactic disks has 
been observed both in external gal axies (e.g., Elmegreen et al.l 
1 19921) and N-body simulations (e.g..lRau tiainen and Said 19991) . 
Given that more than two-thirds of disk galaxies, including our 
own Milky Way, contain both central bars and SS, it is imper- 
ative to establish a strong understanding of the implications of 
this mechanism. Therefore, in this work we study a range of fully 
self-consistent, Tree-SPH simulations, as well as high resolution 
pure N-body simulations, searching for the signature of spiral- 
bar interaction in the distribution of angular momentum: MFIO 
predicted that when both bar and spirals are present in the disk, 
one finds that the changes in angular momentum form a bimodal 
distribution, with two local maxima close to the corotation and 
outer Lindblad resonance of the bar, regardless of the pattern 
speed of the spiral. 



2. Models and results 

While not self-consistent, test-particle simulations allow for a 
full control over the simulation parameters, such as the am- 
plitudes and pattern speeds of bar and SS, and still provide a 
good approximation to self-consistent simulations. Employing 
this method in MFIO, we were able to suppress the effect of tran- 
sient spirals and thus identify the non-linear effect of resonance 
overlap. We showed that the most important signature of this 
mechanism was a bimodality in the changes in angular momen- 
tum in the disk, AL, with maxima near the bar's corotation and 
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Fig. 1. First row: Stellar disk number density contours of the isolated giant Sa galaxy simulation in GaUVIer for 5 time outputs as 
indicated in each panel. Second row: Changes in angular momentum, AL as a function of the initial angular momentum, L. Both 
AL and L are divided by the asymptotic circular velocity Vc, thus L gives approximately the galactic radius in kpc. The locations 
of the bar's corotation and 2:1 OLR are indicated by the dotted and solid lines, respectively. A bimodal distribution indicating the 
work of resonance overlap of bar and spirals is clearly seen in each panel. Third row: The evolution of the radial profiles (left) and 
metallicities (right) for the stellar and gaseous disks. The initial stellar and gaseous disk scale-lengths are indicated by the dashed 
lines. The time steps shown are as in the first row, indicated by solid red, dotted orange, dashed green, dotted-dash blue, and solid 
purple lines, respectively. Fourth row: Distribution of birth radii for stars ending up in the annuU indicated by the green fines (600 
pc) at t=2.5 Gyr 



its outer Lindblad resonance (OLR) regardless of the SS pattern 
speed. 

Hereafter, we analyze self-consistent simulations with strong 
bars and spirals, therefore mixing from both transient spirals 
(SB02) and resonance overlap (MFIO) is expected. However, we 
hope to be able to identify the latter mechanism by detecting 
the aforementioned bimodality, given that the predicted distribu- 
tion of AL generated by transient spirals without a bar is rather 
smooth and not bimodal (see SB02). 



2. 1 . Tree-SPH simulation of a giant Sa spiral galaxy 

We test the predictions of MFIO by analyzing fully self- 
consistent simulations of isolated disk galaxies from the GalMer 



databas jl including a gas component as w ell as star formation 
( Di Matt eo et all l2007l:IChifingarian et alibolOh . 

In Galmer, for each galaxy type, the halo and the op- 
tional bulge are modeled as Plummer spheres, with characteristic 
masses Mh and Mg, and characteristic radii rfj and rg, respec- 
tively. Their densities are given by: 



PH{r) 



and 



PB{r) 



3Mh 



3M« 



1-^ 



-5/2 



-5/2 



(1) 



(2) 



http://galmer.obspm.fr 
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Fig. 2. Same as the second row of Fig. [T]but for all isolated giant GalMer galaxies. Contour levels are the same for all plots and 
the time is t=l Gyr. The locations of the bar's corotation and OLR are indicated by the dotted and solid lines, respectively. The 
rightmost panel shows the only disk lacking a central bar. In this case, resonance overlap from strong multiple spirals gives rise to a 
smooth AL distribution. 



On the other hand, the gaseous and stellar disks follow 
Miyamoto-Nagai density profiles with masses Mg and M» and 
vertical and radial scale lengths given, respectively, by hg and 
Gg, and lit and a*: 



( h/Mg 

4-K 



QgR + (flg + 



at + 



An 



3^z^+hg^)(ag+ ^z^+hg^"^ 



(3) 



(4) 



Here, we first analyze the properties of the isolated giant Sa 
(gSa) galaxy model: it has a total mass of ~ 2.4 x 10''M5„„, 
with Mh = 1.15 X 10"M,„,„ Mb = 0.23 x 10"M,„„, Mg = 
0.09 X 10"M,„„, and M, = 0.92 x 10"M,™„. The scale radii are 
th = 10 kpc, rg -2 kpc, hg - 0.2 kpc, Ug - 5 kpc, /i, = 0.5 kpc, 
and a* = 4 kpc. The initial Toomre parameter of both stars and 
gas is taken to be g = 1 .2 as the initial condition of the Tree- 
SPH simulations. Fig.[T]shows the properties of this isolated gSa 
model as a function of time. 



2.1.1. Identification of tine bimodality 

In the first row of Fig. [1] we plot stellar number density con- 
tours for five time outputs up to 3 Gyr as indicated in each panel. 
We note that the initially strong SS disappears by the end of the 
simulation. The expansion of the disk is already an indication of 
radial migration. 

The second row of Fig. [T] shows contours of the changes in 
the angular momentum of stars, AL, versus the initial angular 
momentum, Lq. Both AL and L are divided by the asymptotic 
circular velocity Vc ~ 280 km/s, thus L gives approximately the 
galactic radius in kpc. The bar's corotation and 2:1 OLR are in- 
dicated by the dotted and solid lines. A bimodal distribution be- 
comes apparent as soon as the bar and SS form (f k 0.2 Gyr). 
While the peak at Lq ^ 4.5 is caused by the bar's corotation, the 
one at Lq ~ 8.5 could oyily be produced by resonance overlap 



(see Fig. 2 in MFIO). At all later times this clear bimodality re- 
mains in the distribution of AL, thus providing strong evidence 
that the effect we see is caused by the spiral-bar interaction. We 
note also that the mixing timescale here is simply too short to be 
the effect of transients only. 

2.1 .2. Extension of tine disk 

The third row of Fig.[T]shows the temporal evolution of the radial 
density (left) and metallicity (right) profiles for both the stellar 
and gaseous disks. The time steps shown are as in the first row, 
indicated by solid red, dotted orange, dashed green, dotted-dash 
blue, and solid purple in increasing order. 

Owing to the vigorous mixing, the disk grows rapidly with 
time displaying a break in the density profile that moves progres- 
sively to larger radii. At the final time (purple lines), the stellar 
and gaseous disks have extended to about ten and seven scale- 
lengths, respectively, while roughly preserving their exponential 
profiles. This is accompanied by a strong flattening in the metal- 
licity gradients, where at t=3 Gyr one can clearly see a reversal 
in the metallicity gradient of stars in the outer region of the disk 
(~ 30 kpc). The stellar metallicities at that location are similar 
to the initial ones at 8 kpc, suggesting that stars in the outskirts 
of the disk originate from regions near the bar's 2: 1 OLR, where 
the changes in AL are most prominent. 

2.1.3. Birtli radii of stars 

The bottom row of Fig.[T]shows the distribution of birth radii for 
stars ending up in the annuli indicated by the green lines (600 pc) 
at f = 2.5 Gyr The black and red histograms show stars on nearly 
circular orbits (of velocity dispersion smaller than 10 km/s) and 
the total population, respectively. In all cases, a large fraction 
of stars are found to originate from the bar's corotation (dot- 
ted Une), as well as an increasingly large fraction from the bar's 
OLR (solid line) as larger final radii are sampled. 

2.2. Tree-SPH simulations of otiier galaxy types 

To assess how our results depend on the peculiar galaxy type 
we considered, we plot in Fig. |2] the changes in angular mo- 
mentum for the stellar disk compo nents of all GalMer iso- 
lated giant galaxy models available dChilingarian et al.l 1201 Ol 
Table 1). Contour levels are the same for all plots and the time 
is t=l Gyr. All apart from the gSd (bulgeless) model develop 
long-lived, central bars at the beginning of the simulation and 
consequently show a bimodality in AL. When we compare the 
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Fig. 3. Power spectrum displaying the pattern speed of the m = 2 
and m = 4 Fourier components, during the entire simulation, 
of the gSa (top) and gSb (bottom) models. Contour levels are 
indicated for each panel. Note that the gSa bar evolves (slows 
down and extends) much more and its two-armed SS is twice 
as strong as for the gSb model (top left). Despite the slightly 
stronger gSb bar, the effect on AL in the gSa's outer disk is much 
more prominent (Fig|2]i because of the combined effect of the bar 
and stronger spirals. 




Fig. 4. Contrasting the effect on AL for a low-resolution (1 .27 x 
10^) and a high-resolution 10^ pure N-body simulations. 



gSb bar, the effect on AL for the gSa is much more prominent 
in the outer disk (FigUJi because of the combined effect of the 
bar and stronger spirals. Strong spirals are thus needed for this 
migration mechanism to be efficient. However, the relation be- 
tween the strength of the perturbers and the angular momentum 
redistribution is nonlinear as previously shown by MFIO. 



gSO and gSa models, which have identical initial conditions ex- 
cept that gSO is gas-free, we find that AL increases by ~ 20% 
when gaseous disk component is introduced. The stability of 
the gas, however, does not appear to play a major role. Only 
a small difference occurs when the gas is made more unstable 
when Qgas is lowered from 1.2 to 0.3 (second and third pan- 
els), where ggns is the Toomre instability parameter for the gas. 
The weaker bar in the gSb model (fourth panel) results in re- 
duced mixing as we predicted in MFIO. Unlike all other models, 
the gSd simulation (rightmost panel) does not form a stable bar. 
The strong mixing obs erved here result s from the large num- 
ber of multiple spirals jELmegreen et all 119921 : iRix and RiekeL 
119931: iMinchev and Ouillenl |2006|) propagating simultaneously 
through the disk during the first ~ 500 Myr. Due to the various 
SS pattern speeds, there are no radii at which the effect of res- 
onance overlap differs distinctly from the rest of the disk (such 
as the bar's corotation and OLR). Consequently, we observe a 
smooth distribution of AL in contrast to the case of bar + SS. 

2.3. Pattern speeds 

Fig.[3]shows a power spectrum displaying the pattern speeds of 
the m - 2 and m - 4 (two-armed and four-armed structure) 
Fourier components, over the whole simulation, of the gSa (top) 
and gSb (bottom) models. Contour levels are indicated for each 
panel in units of Qr, the ratio of the maximum tangential force 
to the azimuthally averaged radial force at a given radius. For 
both simulations, the SS (r ^ 5 kpc) spans a range of pattern 
speeds, always slower than the bar. As expected from the re- 
sults of MFIO, the SS pattern speed has little influence on the 
location of the maxima in the bimodal distribution of AL. Note 
that the gSa bar evolves (slows down and extends over time) 
much more, and that its two-armed SS is twice stronger com- 
pared to the gSb model (top left). Despite the slightly stronger 



2.4. High resolution N-body simuiations 

To determine whether the strong mixing we described above is 
caused by a deficiency in the resolution of the GalMer simu- 
lations, we ran pure N-body collisionless simulations with 10^ 
and 1.27 x 10^ part icles in the disk. Full description and de- 
tails can be found in lWozniak and Michel-Dansac (20()9). With 
a scale-length of 3.5 kpc, the initial conditions of the disk are 
comparable to those of the gSb GalMer model (Fig. |2| but lack 
a gas component. These simulations have no halo, which causes 
a ~ 30% drop in the RC at ~ 10 kpc. 

Fig. |4] shows the changes in angular momentum at t=l and 
3 Gyr for 1.27 x 10^ (top) and 10^ (bottom) disk particles. At 
the beginning of the simulations, the effect is stronger for the 
low-resolution case due to the faster bar formation. However, at 
f Si 2 Gyr both the high- and low-resolution runs yield a simi- 
lar result. This suggests that the magnitude of the radial migra- 
tion induced by the bar-spiral interaction and the timescale over 
which it is effective are not strongly dependent on the numerical 
resolution of the simulations. We conclude that the properties of 
the mixing, i.e., its magnitude and timescale, that we infer from 
the GalMer simulations are only weakly affected by the insuffi- 
cient resolution. 

2.5. Migration in low-mass galaxies 

We now investigate whether the resonance overlap mechanism 
is also efficient in low-mass galaxies. Gogarten et al. (201^ 
showed that although transient spirals can explain extended disks 
for MW-type galaxies, this mixing mechanism is inefficient for 
galaxies with RCs of V,- ~ 100 km/s, such as NGC 300 and 
M33. Nevertheless, both NGC 300 and M33 are observed to 
have extended radial profiles of up to ten scale-lengths. We now 
consider the barred, dwarf Sa (dSa) simulation in GalMer (RC 
of 100 km/s and initial scale-length of 1.3 kpc) to see how its 
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Fig. 5. The time evolution of the stellar disk density (left) and 
metallicity (right) profiles for the dwarf Sa isolated GalMer 
model. Time outputs and colors are the same as the third row 
of Fig. [1] This is a low-mass galaxy with a rotation curve of 
~ 100 km/s similar to NGC 300 and M33. The initial disk scale- 
length is indicated by the dashed line. The disk extends to over 
ten scale-lengths in less than 3 Gyr. 

density distribution and metallicity gradient evolve with time. 
In Fig. |5] we plot the time development of the stellar disk den- 
sity and metallicity profiles, for the same time steps as in Fig.[T] 
We can clearly see that the stellar disk extends to more than ten 
scale-lengths in ~ 3 Gyr, while preserving its exponential mass 
density profile. At the same time, the metallicity gradient be- 
comes flat in less than 1 Gyr. Both the predicted extent of the 
stellar disk and the flattening of the metallicity profile are con- 
sistent w ith what has been recently reported for nearby low -mass 
galaxies dBland-Hawthorn et al. l l2005HVlaiic et al.Ll2009h . 

3. Conclusions 

We have examined the redistribution of angular momentum in 
galactic disks by means of Tree-SPH and high-resolution pure 
N-body simulations. We have found that resonance overlap of 
multiple patterns (such as bar -i- SS or SS + SS) induces strong 
exchange of angular momentum throughout the disks in agree- 
ment with the predictions of MFIO. Since in the self-consistent 
simulations analyzed in this work spirals may be transient, we 
should also expect a contribution from the SB02 radial migra- 
tion mechanism. However, for the short timescales considered 
here (< 1 Gyr) transients would simply have a very small ef- 
fect (SB02). The resonance overlap mechanism is clearly iden- 
tified by a bimodality in the changes in angular momentum, AL, 
caused by the bar's corotation and 2: 1 OLR (Fig.[Tl|2] andHJi. We 
have contrasted this to a simulation lacking a stable central bar 
(Fig. |2] rightmost panel), where the AL distribution is smooth. 
The effect is especially strong when a gaseous component is 
present, as a result of the exchange of L between gas and stars 
coming from the gravity torques, the ph ase shift between the two 
components, and the gas dissipation (iBournaud and Combesl 
l2002h . Depending on the amount of gas and the strength of the 
bars and spirals, the metallicity gradients can flatten in less than 
1 Gyr (Fig. [U. This is in drastic contrast to the current under- 
standing that galactic disks need a Hubble time for sufficient 
mixing (SB02. lRoskar et al.ll2008h . 

How can we tie our re sults to galactic disk evolution? 
iBournaud and Combesl (l2002h followed the detailed processes of 
bar formation, bar destruction, and bar re-formation, while vary- 
ing the disk to bulge ratio. These authors identified three bar for- 
mation episodes in a Hubble time. In this model, we can regard a 
given GalMer isolated galaxy simulation as one such episode of 
a gas accretion event during a galaxy lifetime. In this scenario, 
the rapid flattening of the metallicity gradients expected from the 
vigorous migration (Fig.lTJ would be followed by a gas enrich- 
ment at each bar re-formation event resulting in the rebuilding of 



the gradient. The presence or not of a metallicity gradient, or its 
intensity, would then be an indicator of the bar/accretion phase 
of the galaxy. 

The mechanism described in this study works even in low- 
mass galaxies (Fig.|5]) and can thus for the first time provide an 
explanation for the e xtended disk profiles observed in galaxies 
with Vc ~ 100 km/s dBland-Hawthorn et al.L I2005L Minchev et 
al. 201 1, in preparation). 
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